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ABSTRACT 


The presence of flare related acoustic emission (sunquakes) in some flares, 
and only in specific locations within the flaring environment, represents a severe 
challenge to our current understanding of flare energy transport processes. In 
an attempt to contribute to understanding the origins of sunquakes we present a 
comparison of new spectral observations from Hinode’s EUV imaging Spectrom¬ 
eter (EIS) and the Interface Region Imaging Spectrograph (IRIS) of the chro¬ 
mosphere, transition region and corona above a sunquake, and compare them to 
the spectra observed in a part of the flaring region with no acoustic signature. 
Evidence for the sunquake is determined using both time-distance and acoustic 
holography methods, and we find that, unlike many previous sunquake detec¬ 
tions, the signal is rather dispersed, but that the time-distance and 6 and 7 mHz 
sources converge at the same spatial location. We also see some evidence for 
different evolution at different frequencies, with an earlier peak at 7 mHz than 
at 6 mHz. Using EIS and IRIS spectroscopic measurements we find that in this 
location at the time of the 7 mHz peak the spectral emission is significantly more 
intense, shows larger velocity shifts and substantially broader prohles than in the 
location with no sunquake, and that there is a good correlation between blue- 
shifted, hot coronal, hard X-ray (HXR) and red-shifted chromospheric emission, 
consistent with the idea of a strong downward motion driven by rapid heating 
by non-thermal electrons and the formation of chromospheric shocks. Exploiting 
the diagnostic potential of the Mg H triplet lines, we also find evidence for a 
single, large temperature increase deep in the atmosphere, consistent with this 
scenario. The time of the 6 mHz and time-distance peak signal coincides with 
a secondary peak in the energy release process, but in this case we hnd no ev¬ 
idence of HXR emission in the quake location, but very broad spectral lines. 



strongly shifted to the red, indicating the possible presence of a signihcant flux 
of downward propagating Alfven waves. 


Subject headings: 
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Introduction 


In the ’standard’ 2D eruptive flare paradigm (for a review see IShibata fc Magara 


(1201 ih ) energy is released via magnetic reconnection in the solar corona where it is converted 


to heat, mass motions, and the kinetic energy of accelerated particles, and is subsequently 
transported both out into the heliosphere and down into the deeper solar atmosphere. 
Observations of the subsequent response of the plasma to this energy release do, in some 
aspects (e.g the appearance of flare loops), appear consistent with expectations based on 
this paradigm, but it is clearly unable to account fully for all the observed elements of the 
flare process. 

One area that represents a significant challenge for the ’standard’ model is the 


observation of localized acoustic dis 


sunquakes (jKosovichev fc Zharkova 


urban ces in response to flaring activity. These so-called 


19981 ) arguably represent one of the most severe 


ar atmosphere, despite the rel atively smal 


constraints for energy transport through the so 
fractio n of total flare energy involved (< 10“^, iDonea &: LindsevI (1200511 ). As 


Judge et al. 


(1201411 point out, if we assume (as the standard model tells us) that the flare energy 
is released in the corona, then in order to drive an acoustic disturbance in the solar 
interior, the energy must propagate through nine pressure scale heights. Since acoustic 
enhancements are also not detected for every flare, and only in certain locations within the 
flaring active region when they do occur, understanding their origins can also provide new 
insights into the effects of environmental (e.g. plasma, magnetic held) conditions on flare 
energy transport. 


The first sunquake was observed by 


Kosovichev fc Zharkoval (Il998ll . and most 


subsequent studies have indicated that the sunquake originates in the impulsive phase of 


the flare and is gener ally well-correlated with t 


l e locations of enhanced white-light and 


(Zharkova & Zharkov 

20Q7; 

Donea & Lindsev 

LO 

O 

O 

CM 


200511. This observed 
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correlation has led to suggestions that chromospheric shocks, which arise as the result of 
the pressure transients driven by the hydrodynamic response of the ambient plasma to 
the precipitation of energetic particles into the chromosphere; pressure transients that are 
related to back-war ming of the photosphere by enhanced chromospheric radiation, (e.g. 


Donea fc Lindsey! (120051 )): and Lorentz force tr ansients that occur whilst the coronal held 


is being reconhgured (e.g. 


Hudson et ah 


f)2012l) ). may be potential drivers of the acoustic 


emission. However, th e launch of SDO, coupled with now hnely-honed detection techniques 
flZharkov et al.ll2011bh . has led to a series of sunquake observations that show signihcantly 
weaker correlations between the location of the acoustic sources and the areas of strongest 


connection to the ends of the erupting hux rope (I 

^edram & Matthews 

2012 

2nna: 

Matthews et al. 

2011; 

Zharkov et ahl 

201 3al 

b) suggesting that we are 


Zharkov et ah 


from understanding the mechanism ( s) resp onsible for producing these important events. 


Indeed, recent work by 


Judge et al. 


(120141) reports unique and important observations in 


the infra-red (IR) during a sunquake that constrains the depth of the hare heating to be in 
the photosphere (lOOilOO km), but also demonstrates that the level of non-magnetic and 
magnetic energy huxes consistent with their observatio ns are too low to dr ive the acoustic 


transient. In this work we examine the same hare that 


Judge et ahl (120141) analysed with 


a focus on complementary spectroscopic observations of the corona, transition region and 


(Culhane et al. 

2007 

) and IRIS ( 

De Pontieu et al. 

2014) 


20141) in the location above the quake 


with spectra above a region in the Northern hare ribbon that has no acoustic enhancement 
in order to look for signatures of downward propagating energy in these regions and any 
diherences that may shed additional light on the sunquake origin. 
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2. Observations and Data Reduction 


On 29 March 2014 an XI flare (SOL2014-03-29T17:48) occurred in active region NOAA 
12017. This was an exceptionally well-covered event, observed from space by Hinode; the 
Solar Dynamics Observatory (SD O); IRIS, and the Reuven Ramaty High Energy Solar 


Spectroscopic Imager (RHESSI) flLin et ahll2002h . as well as from the ground by t 


he FIRS 


(Faci lity Infrared Spectrometer) and Imaging Bldimensional Spectrometer (IBIS) fjCavallini 


20061 1 at the Dunn Solar Telescope (DST). The photospheric sig natures of th i s even t, and 


their relation to the sunquake, have been extensively studied by 


Judge et ah 


(1201411 . and 


the work we present here is complementary in that it probes conditions in the upper layers 
of the atmosphere. Thus, we present observations primarily from Hino de EIS, IRIS, SDO’ s 


HMI and RHESSI, with SDO’s AIA (Atmospheric Imaging Assembly) ( Lemen et al.l 12012 1 


providing supporting context. The temporal evolution of the flare as observed by GOES and 
RHESSI is plotted in the right-hand column of Figure 1, with the vertical lines indicating 
the time intervals of the sunquake, while on the left is the IRIS Mg H slit-jaw image from 
the peak of the impulsive phase, illustrating the ribbon structure in the chromosphere. 
Marked on the slit-jaw image are the location of the IRIS spectrograph slit (vertical line), 
and the locations of the sunquake (SQ), and of a region within the northern ribbon that 
we call the non-quake (NQ) region. This region was chosen for comparison because it 
represented an area of very strong UV and EUV intensity enhancement outside the SQ 
location, and was scanned by the IRIS spectrograph. 


EIS observed the region from 14:05 - 17:57 UT, executing the high cadence flare study 

n n n n n 

PRY_flare_l, which covered an FOV of 44 x 120 centred on (499 ,274 ) with the 2 slit in 

n 

4 steps to produce a raster cadence of 134 seconds. The study included eight line windows, 
including the He H 256.3 A (log Te=4.7)and Fe XXHI 263.76 A (log Te=7.1) lines, which 
are the main focus of this work. The EIS data were calibrated using the standard software 
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routine eis_prep, and the slit tilt and the orbital variation of the line position were corrected 
for. The lines were htted with either a single or two-component Gaussian £t, which allowed 
us to construct intensity, Doppler velocity and width information for each raster. 

SDO’s HMI and AIA instruments both observed the whole disk of the Sun and the flare 

throughout its duration. Standard calibration (aia_prep) and image alignment procedures 

were applied, and the other data were aligned to these data. Co-alignment with EIS can 

be challenging given the small FOV, and the temporal smearing that is introduced by the 

rastering of the slit. Thus, while the offset between EIS and AIA is now routinely recorded, 
// 

and is better than 5 in x and y, we confirmed our alignment using feature matching of 

the integrated intensity in the Fe XII 192.41 A window and images from the AIA 193 A 

bandpass recorded at the mid-point of the EIS raster. This alignment, which we estimate 
n „ 

to be good to ± 2 , then allowed us to align EIS, AIA 1600 A and the IRIS slit-jaw data. 


IRIS observed the event both with the slit-jaw imager and the spectrograph. The 

// n n // 

spectrograph covered an FOV of 14 x 174 centred at (491 ,282 ) with 8 slit positions at a 
9.4 second exposure time, giving a total raster cadence of 72 seconds. The line list included 
the Mg II h & k lines at 2803 and 2796 A respectively, the Si IV 1393.76 and 1402.77 A lines 
and the C II 1334.53 and 1335.71 A lines. The slit-jaw imager provided images in the 1400 A 
and 2796 A bandpasses covering a FOV of 167 x 174 from 14:09 -17:54 UT. The IRIS pixel 


size for both the slit-jaw imager and spectrograph is 0.166 . We primarily used the IRIS 
level 2 spectrograph data in our analysis, and the level 2 calibration includes corrections for 
dark current, flat-fielding as well as geometric and wavelength corrections. Although IRIS 
observed in several wavelength windows, the intensity of the flare was such that many of 
the spectra were saturated during the peak of the impulsive phase. This inher ently limits 


our an alysis and interpretation in a way that the ground-based observations of 


Judge et ah 


(120141) are not subject to, and for this reason we note the importance of combining both 




ground and space-based data sources when possible. In this work we focus primarily on the 
Mg II h and k lines, and the Si IV line. Both the Mg II and Si IV do show some saturation 
during the flare, and we exclude these data from our spectral htting. However, since these 
wavelength regions are inaccessible from the ground, these data still provide important 
information about the flare process, and despite the saturation, substantial changes in the 
Mg II h and k lines are seen during the flare, including large velocity shifts, line broadening, 
disappearance of the central reversal and the appearance of the subordinate 3p^P -3d^D 


triplet lines, which offer diagnostic potent ial for flare heating in the l ower atmosphere t 


lat 


has so far been exploited relatively little flFeldmau fc Doschek 


19771: 


Pereira et ah 


2ni5h 


Both the IRIS and EIS spectrographs were operating in a scanning mode during the event 
studied, and the time to build a spectral image using rastering is inherently longer than the 
time taken to produce an image using broad-band hlters. However, as can be seen from 
hgures 4-8 and 10 - 12, we have still been able to sample spectra at important points 
during the flare evolution. 


RHESSI observed for the entire duration of the flare and the data were used to 
construct light-curves in a range of energy bins, and to reconstruct images in the 25-50 and 
50-100 keV energy ranges. For the image reconstruction we used the PIXON algorithm 
with a 20 second integration time. Figure 2 shows the location of the reconstructed HXR 
sources at 25 -50 (purple) and 50 - 100 keV (black) with respect to flare ribbons observed 
in the IRIS Mg H slit-jaw images. 


2.1. Acoustic methods 


by 


As noted in the i ntroduction, the presence of the sunquake has already been reported 


Judge et al. 


fl2014l) . and the holography data that we present here are independent 


conhrmation of those results. However, we also include complementary time-distance 
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RIS SJI SDO 29-Mar-2014 17:46:04.710 UT 



460 480 500 520 540 560 580 

X (arcsecs) 



17:30 17:45 18:00 18:15 18:30 

Start Time (29-Mar-14 17:29:57) 



17:30 17:40 17:50 18:00 18:10 18:20 18:30 
Start Time (29-Mar-14 17:30:00) 


Fig. 1.— Left: IRIS Mg II slit-jaw image showing the location of the sunquake (SQ) and 
the non-quake (NQ) region and the location of the spectrograph slit (vertical line). Right: 
top: GOES X-ray intensity for the flare; bottom: RHESSI count-rates in the 12-25 keV 
(red) and 25-50 keV (blue) energy ranges. The solid back (dashed orange) lines represent 
the time intervals of strong photospheric Doppler velocities and the derived sunquake onset 
time, respectively. 
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analysis that confirms the location and timing of the snnqnake. For the holography 


we also nsed t 


flScherrer et al 


r e HM I photospheric data obtained aboard SDO with 45 second cadence 


20121) . For helioseismic analysis the data are tracked at the Snodgrass 


differential rotation rate and re-mapped nsing Postel projection for helioseismic purposes. 
In order to detect acoustic emission associated with the flare, we applied both acoustic 
holography and surf ace time-distance a nalysis to the tracked Dopplergram data. The 


holography method flDonea et al. 


19991) uses Green’s function G+dr — r'|,f — t') which 


prescribes the acoustic wave propagation from a point source, to essentially ’backtrack’ the 
observed surface signal This allowed egression images to be reconstructed showing 

the sub-surface acoustic sources and sinks: 


/ dt' / d'^r'G+{\r — r'\,t — ,t'), 

J Ja<\r—r'\<b 


( 1 ) 


where a, b dehne the holographic pupil. The egression power is then dehned as : 

P{z,r) = j dt\Hj^{r,z,t)\^dt (2) 

where r' and 6 are the polar coordinates describing the vector |r — r'| at the solar surface. 


The time-distance analysis flKosovichev fc Zharkova 


1998 


Zharkov et al. 


2011b 


2013b|) uses the surface line-of-sight velocity obtained from the three-dimensional series 
of observations (two horizontal dimensions plus time) by selecting the source location, 
re-writing the observed surface velocity signal in polar coordinates relative to the source, 
then contracting the azimuthal dimension via integration applying a Fourier transform. 
The result is a time-distance diagram presenting the radial surface distance from the source 
versus time seen as a ridge following the time-distance relationship for the acoustic waves 
derived from Dopplergrams, pointing to the circular ripples propagating on the solar surface 
from the source where a flare occurs. 

Photospheric velocity changes are seen in the HMI Doppler data within the southern 
flare ribbon at (519 ,262 ) between 17:45:31 - 17:47:46 UT. During large flares it is likely 
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that HMI channels will be affected by intensity core reversals, a nd indeed t 
the magnetic held data. Also, the photospheric data analysed bv ljudee et al 


lese are seen in 


(2014) indicate 


qnite modest velocities in the quake location. Most of the sunquake detections made using 
acoustic holography reported up until recently have been characterised by a well-dehned 
localised acoustic kernel. Such kernels are normally seen in egression measurements at more 
than one frequency. While the shape and extent of the kernels may vary from frequency to 
frequency, they tend to be cospatial. The shape, longevity and prominence of kernels may 


also depend on the choice of the pupil used in egression. B ut again, the 


emission remains approximately the same. As Figure 1 in 


Judge et al. 


ocati on of strongest 




fl2014l) shows, and is 


independently verihed by our measurements, the 6mHz egression data in this event are quite 
noisy, indicating several potential kernels. We also found that the shapes and locations of 
these kernels varied considerably from frequency to frequency as well as with the choice of 
pupil. To localise the egression power emission, we used a 15 by 15 Mm integration box to 
analyse the acoustic emission statistics for all data points at different frequencies, marking 
the locations where signal passes 5(7 threshold. This, combined with time-distance analysis, 
allowed us to verify and localise the source. In Figure 3, we show the 6 and 7mHz egression 
power contours (cyan and magenta, respectively), together with the 25 -50 and 50 -100 
keV HXR contours overlaid on the IRIS 2976 A slitjaw image at 17:46 UT, indicating that 
there was strong energy deposition in this location at many wavelengths. Figure 2 shows 
the time-distance plot associated with the acoustic source at (519 ,262 ), and from this 


we determine an onset time is in t 
is also the focus of the IJndge et al 


le range 17:48:31 - 17:51:31 ±45 seconds. This source 


(I 2 OI 4 J ) work, and they report two onset times for the 


egression sources at 5.5 and 6 mHz of 17:46 and 17:51 UT, respectively. In our analysis we 
also hnd some evidence for different evolution at different frequencies, with a peak time 
between 17:47:46 - 17:50:01 UT at 7 mHz, and between 17:50:46 - 17:51:31 UT at 6 mHz. 
It should be noted that because of the hltering involved, the uncertainties in onset times 
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derived from holography can be ±5 minutes. 


3. Results 


The egression source located at (519 ,262 ) has the greatest statistical signihcance, 


an d is also the 
of I, Judge et ah 


ocati on of the time-distance source. In order to complement to the work 


( 20141 ) in trying to understand the processes that led to the production 


of the acoustic disturbances at this particular location of the flaring environment, we use 
spectroscopy to probe the conditions that prevailed in the overlying corona, transition 
region and chromosphere at this position, and compare them to conditions at another 


bright location (NQ) in the northern flare ribbon at (511 ,272 ). These positions can be 
seen in Figure 3 to correspond to the sites of strong Mg II emission during the peak of the 
impulsive phase at approximately 17:46:30 UT, and for the sunquake, to the strongest 6.0 
mHz egression contours and HXR emission. The NQ region is at the edge of the strongest 
HXR emission in the northern ribbon, and this location is at the Eastern edge of the FOV 
of the IRIS spectrograph, so we were unfortunately unable to study the Mg II spectral 
response any closer to the centroid of the HXR emission in the Northern ribbon. 


The intensity and dynamics of the flare produced complex prohles in all the spectral 
lines that we examined, as can be seen in hgures 4, 5, 6 and 7, where we plot the He H 
256.3 A Fe XXHI 263.76 A Si IV 1402.8 A and Mg H h and k spectra for the NQ location 
(left column) and SQ location (right panels), respectively. 


3.1. He II spectra 

The He H 256.3 A spectra were htted with a combination of single and double Gaussian 
component hts, and the reference wavelength was derived from an average of the pre-flare 
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spectra to be 256.36 Ain Figures 4 and 5 we plot the He II spectra from the NQ and SQ 
locations (c.f. Figure 1 for locations within the flare ribbons). Figure 4 shows spectra from 
the impulsive phases of the flare (14:42 - 17:47 UT), while Figure 5 shows the spectra from 
the gradual phase (17:48 - 17:55 UT). Since both EIS and IRIS raster in order to build up 
an image, the timing at each slit location within the image differs by the time taken to move 
between raster positions; in this case 24 seconds, and so it is important to keep in mind 
that each location is not measured simultaneously. However, it can be seen from this hgure 
that the He H prohle in the SQ location is consistently broader and of greater intensity 
than in the NQ location, and that it is red-shifted with respect to the reference wavelength, 
with the greatest red-shift of ~80 kms“^ occurring at 17:46:38 UT. It is interesting to note 
that an extended red ’tail’ begins to develop in the SQ location around 17:51 UT, with a 
very strong and highly shifted red component apparent at 17:55 UT. The spectra in the NQ 
location indicate a signihcant red-shift and red wing component at 17:47 UT, but thereafter 
the line centroids indicate only small deviations from the reference wavelength, indicating 
that there are no signihcant motions of the emitting plasma. 


3.2. Fe XXIII spectra 

Prior to 17:46 UT, there is no signihcant Fe XXHI 263.78 A emission visible in 
either the SQ or NQ locations, although it is observed further East in the active region. 
However, as can be seen from Figure 6, the spectra in both locations show a substantial 
blue component from 17:46 UT onwards, reaching ~280 kms“^ in the SQ location. In 
the NQ location, low intensity blue wing emission is seen until 17:51 UT, while in the SQ 
location the prohle is symmetric but shifted to longer wavelengths from 17:48 UT onwards. 
The peak red-shift velocity in the SQ location is 26 kms“^ at 17:51 UT. 



14 


3.3. Si IV spectra 


From pre-flare spectra we derive a reference wavelength of 1402.82 A for the Si IV 
line. As can be seen from Figures 7 and 8, where we plot the NQ and SQ spectra together 
with reference light-curves from GOES and RHESSI, it is clear that both the NQ and SQ 
locations suffer from saturation at a number of times during the flare impulsive phase (the 
third panel down on the right in figure 7 is blank for this reason). However, the NQ location 
shows strongly blue-shifted compo nents at 17:43 and 17:44 UT, probably associated with 


the onset of the hlament eruption flKleint et ah 


20151) . followed by saturation during the 


flare peak (17:45-17:46 UT), and the subsequent development of a red ’tail’ from 17:48 
UT onwards. In the SQ location, there is an initial red-shift at 17:43 UT, which develops 
into an extended red tail at 17:44 UT. The spectra are saturated at 17:46 and 17:47 UT, 
and we then see the development of a strong red component from 17:48 UT, with a peak 
downward velocity of 86 kms“^ at 17:52 UT. In Figure 9, we plot the intensity and Doppler 
widths derived from the Si IV line hts, omitting the saturated spectra. While the intensity 
and Doppler width at the flare peak are not accurately reflected in these plots, what is 
clear is that there are multiple peaks in the SQ location, with the first occurring before the 
impulsive phase peak, at ~17:43 UT, and a later strong burst after the impulsive peak at 
~17:51 - 17:53 UT. A substantial increase in Doppler width is also seen coincident with the 
later intensity increase. 


3.4. Mg II h and k spectra 

As for the Si IV line, the Mg II spectra are heavily saturated in the location identified 
with the sunquake, as can be seen in the third panel on the right in Figure 10. This made it 
impossible to fit the spectra and derive the intensity or velocity at this position. However, 
we were able to £t the spectra in the SQ location in preceding and successive rasters (Figure 
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11), as well as the spectra in the NQ location. It should be noted that the Mg II and Si 
IV lines can be optically thick, and hence hts to the spectra do not simply reflect plasma 
properties, but do provide a means to comment on relative changes in the line centroid and 
prohle. The spectra in both locations showed complex variations in the line prohles of both 
the h and k lines, as well as in the prohles of the companion triplet lines that arise from the 
3p^P and 3d^D transitions at 2791.60 A and 2798.75 and 2798.82 A . In particular we note 
the disappearance in the central reversal, the increase in line width that can be observed 
in both the NQ and SQ locations at the time of the hare peak, the increase in intensity 
and width of the triplet lines, and the subsequent narrowing and decrease in intensity of 
both the resonance and triplet lines in the NQ location following the hare peak. In the 
SQ location similar behaviour is observed until 17:51 UT, when there is evidence of the 
development of a red component in the h and k lines, and a decrease in the intensity of the 
triplet lines. The subsequent rasters at 17:52 and 17:53 UT then show very broad h and k 
prohles, and much reduced peak intensity in the triplet lines, although these are also broad. 
The ’jagged’ nature of the h and k lines is puzzling, and may be due to an instrumental 
ehect. However, the observed intensities at this time are well below saturation levels. 

Using the hts to the Mg II h and k lines, we derived the k/h intensity ratio at both 
the SQ and NQ locations from 17:30 - 17:53 UT, as shown in Figure 12. From this hgure 
we can see that the k/h ratio in the SQ location (left column, second panel) shows little 
variation prior to hare onset, maintaining a value of ~1.1, but then begins to rise at hare 
onset (17:44 UT), peaking around 17:50 UT and then falling to pre-hare values. Note that 
because of the saturation this plot does not include an estimate of the ratio at the hare 
peak at 17:46 UT, nor does the top panel include an estimate of the intensity at this time. 
In the NQ location (right column), the pre-hare ratio is similar, with a maximum of 1.37 
observed at hare onset (17:44 UT), followed by a subsequent drop to a minimum at 17:45:36 
UT, and then a gradual return to pre-hare levels. In the optically thin limit one would 
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expect this ratio to be 2, and the ratio at all times prior to and dnring the flare remains 
less than 2, indicating that the lines remain optically thick. The presence of the saturation 
makes it impossible to comment on the behaviour of the ratio at the peak of the impulsive 
phase in the SQ region, but it is intriguing that the ratio decreases significantly compared 
to pre- and post-flare values during the period 17:44 -17:49 UT in the NQ region, but does 
not appear to do so in the SQ region. 


We also measured changes in the Doppler widths at the SQ and NQ locations in both 
the h and k lines, and these are also plotted in Figure 12; again, SQ in the left-hand column 
and NQ in the right-hand column. In this case we have used the extent of the saturation 
at flare maximum to make an estimate of the line width at this time, but this should only 
be viewed as an upper limit. For both locations it can be seen that there is little variation 
in line width prior to flare onset, and the widths are very similar in both locations. A 
substantial increase in line width is seen at flare onset and peak in both cases, but even 
exercising due caution because of the saturation, it is clear from the spectra that the h and 
k line widths at the SQ location are at least 50% larger those seen in the NQ location, 
reaching ~1.6 AAs can also seen in the spectra in figure 7, the h and k line widths then 
increase again in the SQ location at ~17:50 UT. The saturation in the h line is less than 
for the k line, so we use this line to derive an upper limit of 1.833 A for the line width 
at 17:46:04 UT. Subtracting the instrumental and thermal widths, this corresponds to a 
non-thermal velocity of 190.1 k ms~^. If we assume that this broadening represents a flux 
of Alfven waves, then following 


Judge et al. 


(120141) we can estimate the associated wave 


energy density: 


Kave = pVA{Vnt) erg cm ^ s ^ 


(3) 


where p is the mass density, va is the Alfven speed, and Vnt is the velocity amplitude 
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of the Alfven wave. 


Pereira et al 


(120151) find that the Mg II triplet lines are sensitive to 
column densities in excess of 5 xl0“^ g cm~^ , so since we ar e unab le to directly measure 


the density, p, we use the value tabulated by 


Vernazza et al. 


fjl98ll) for a column mass of 


7.82x10 g cm which corresponds to a m ass density of 5.49 xlO ^ g cm Using B 


Judge et a. 


= 800 G, as measured by 
6.04x10^^ erg s“^ cm“^. Ijudee et al 


(1201 


), we calculate the energy flux density to be 


(1201411 calculate a lower limit for the peak of the 


acoustic power to be 5x10® erg s“^ cm“^. We stress that since the spectra are saturated, 
and this may artihcially enhance the line width, the energy flux derived at this time is an 
upper limit, and we also note the large (~140 G )uncertainty on the magnetic held strength 
measurement. However, the spectra at 17:52 UT are not subject to saturation, and a similar 
calculation yields a non-thermal velocity of 152.2 kms“^ and associated wave energy hux of 
3.87 xlO^^ erg s“^ cm“^ at this time. The error on the magnetic held strength measurement 
is not large enough to reduce either of these values by two orders of magnitude. 


The saturation in the resonance lines in the SQ location at the hare peak (17:46 UT) 
makes it impossible for us to derive a Doppler velocity at that time, although it is apparent 
from Figures 10 and 11 that at 17:46, 17:52 and 17:53 UT these lines are either red-shifted 
or have a substantial red wing component. In fact, hts to the k line at 17:52 UT indicate a 
red component with a velocity of > 100 kms”^. The 2798.75 and 2798.82 A triplet pair are 
also saturated, at 17:46 UT, but the line at 2791.60 A is free from saturation, and so we are 
able to estimate velocity for this line. A substantial red component is seen that peaks at 54 
kms“^, with a width of 0.899 A that causes it to blend with the Fe I line at 2792.394 A , 
which also displays red-shift of 54 kms“^. 


In Figure 13 we plot the relationship between the emission at 17:46 UT (left) and 17:52 
UT (right) from different temperature plasmas with respect to the egression power at 6 and 
7 mHz derived from acoustic holography. The background image in both cases (reverse red 
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colour table) is the emission in the blue wing of Fe XXIII at 280 kms“^. On this background 
image we then plot the 6 and 7 mHz egression contours in cyan/magenta for both times; 
the 25-50 and 50-100 keV HXR emission contours at 80 and 90% of the maximum (purple 
and black, respectively) at 17:46 UT, and the 25-50 keV (purple) contours at 17:52 UT; 
emission from He H at 60 and 80 kms“^ (magenta) at both times, and emission from the 
red wing of the Mg H 2791.6 A triplet line at 70 kms“^ at both times. 


Recent work by 


Pereira et al 


fl2015l) highlights the diagnostic potential of the Mg H 


triplet lines, hnding that the presence of these lines in emission requires steep temperature 
increases in the lower chromosphere (> 1500 K) and typical electron densities of 10^® m“® 
or greater. They also find that the ratio of the core to wing intensity is correlated with 
the temperature increase required to drive the lines into emission, and that the shape of 
the emission line profile provides further information about the temperature profile. Their 
simulations indicate that multiple peaked prohles are generally the result of a temperature 
prohle that varies rapidly with height, while lines with a single peak are consistent with a 
single, strong temperature increase at higher altitudes. When there is a strong temperature 
increase low in the atmosphere, and at high density, then their simulations predict that the 
emission will be seen in the far wings of the line, but that the line core w ill show a central 


reversal, and the resonance lines will be very broad. 


Pereira et al. 


(1201511 demonstrate the 


correlation between core to wing intensity and temperature increase for the long wavelength 
pair of the triplet lines. For the 29 March 2014 flare, this line pair is saturated at the flare 
peak. However, we can derive the 2791.6 A core to win g ratio, and w e find this to be 


Pereira et al. 


(120151) . indicates a 


2.21 for the SQ location, which comparing to Figure 5 of 
temperature increase T(Smaa:) - T^m — 2000 K, where T(Smax) is the temperature at the 
height where the source function maximum occurs, and T^m is the temperature minimum 
between column masses of 10 - 0.1 g cm“^. The prohle of the triplet line at 2791.6 A shows 
an asymmetric central reversal, with a strong enhancement in the red wing, while the h and 
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k lines are clearly strongly broadened, suggesting that the temperature increase in this case 


occurring deep in the atmosphere, and at high densities. An alternative exp! 
the increased intensity o f the M g II triplet lines may be found in the work of 


anation for 


and 


Judgel (1200511 


Laming fc FeldmanI (1199211 who discuss the effects of intense shorts bursts of heating 


that create ’ionising’ plasma conditions, whereby there is a strong burst of emission in the 
metastable transition as the result of more effective population of these levels. 


4. Discussion and Conclusions 


We report new spectral observations of the chromosphere, transition region and corona 
from Hinode EIS and IRIS above the location of the sunquake detected during the XI 
flare of 29 March 2014, which complement the observations and analysis already done by 


Judge et al. 


(1201411 on this event. We compare the spectra from the upper atmosphere 


above the sunquake with spectra from another location in the flaring ribbons that exhibits 
strong emission at multiple wavelengths, but which has no associated acoustic disturbance. 
Using holography and time-distance methods, we c onhrm the existenc e of an acoustic 


impact in the photosphere and below, as reported by 


Judge et al. 


( 20141 ). The time-distance 


analysis indicates an onset time between 17:48:31 - 17:51:31 UT ±45 seconds, while the 
holography suggests some differences in evolution with frequency, with peak times between 
17:47:46 - 17:50:01 UT at 7 mHz, and between 1 7:50:4 6 - 17:51:31 UT at 6 mHz, which are 
also consistent with the hndings of 


Judge et al. 


(1201411 . During the time of the 7 mHz peak 


signal we hnd that in the location of the quake there are spatially correlated upflows (~280 
kms“^) in the hot Fe XXHI line, intense HXR emission in the 50 - 100 keV energy range, 
downflows of 80 kms“^ in the He H line, strongly broadened emission in the Mg H h and 
k lines (~1.78 A ) and intense and broad emission in the Mg H triplet lines, including a 
red wing component with velocities in the range ~ 54 - 70 kms“^. In contrast, in the NQ 
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location we find that, while the Mg II h and k and triplet lines are intense and broadened 
with respect to pre-flare values, the intensity increase is an order of magnitude less, and 
the increase in line width is around a factor of two smaller. At the time of the 6 mHz peak 
signal (~ 17:51 UT), while the GOES and integrated HXR data (see Fig. 1) indicate a 
likely second burst of energy release, we hnd no increase in Mg II or Si IV intensity, line 
width or Doppler velocity in the NQ location, but at the SQ location we hnd an increase 
in Mg II and Si IV intensity, strongly broadened Si IV and Mg II h and k prohles (~1.5 A 
), and a red-shifted component of > 100 kms“^ in the Mg II lines. However, in this time 
interval while we observe the Mg H triplet lines to be broadened, the intensity is two orders 
of magnitude less intense than at 17:46 UT, and we observe no measurable HXR emission 
at 25 - 50 or 50 - 100 keV in either the SQ or NQ location. 


We computed the intensity ratio of the Mg H k to h lines in both locations, omitting the 
saturated spectra in the quake location at 17:46 UT. We hnd this ratio to be signihcantly 
less than 2 for both locations, before and during the hare, but with a signihcant impulsive 
increase at hare onset, followed by a sharp decrease. In the SQ location the ratio is 
comparable to pre-hare values at the time of the 6 mHz peak, while it is signihcantly 
lower in the NQ location. Comp aring our observations to the simulations and diagnostics 


presented in 


Pereira et ah 


( 2015 1 we hnd the that prohle of the Mg H triplet line at 2791.6 


A and the ratio of the wing/core intensity in this line are consistent with a single, strong 
increase in temperature of ~2000 K deep in the atmosphere in the location of the quake. 
However, we note that the increase in intensi t y of t hese lines is also consistent with the 


ionising plasma condition described by 


Judgel ( 20051) 


WolhI (jl972[l hrst suggested that solar hares may drive ac oustic disturbances in th e 


solar interior. The subsequent detection of the hrst quake by 


Kosovichev fc Zharkova 


(119981) not only conhrmed these predictions, but suggested that the driver of the acoustic 
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disturbance was a strong compression shock generated as the result of sudden and rapid 


heating of the chromosphere w ithin the framework of the ’thick-target’ model e.g. flBrown 


1971 


Kostiuk fc Pikelner 


19751) . In this scenario, the rapid heating and compression drive a 


radiative instability that forces the heated chromospheric plasma upwards (chromospheric 
evaporation), while generating an opposite downward propagating shock, which then drives 
the acoustic disturbance. In this picture, the heating and subsequent evaporation are driven 
by the non-thermal electrons that are generated by the reconnection process in the corona. 

The majority of sunquakes discovered prior to the launch of SDO showed a good 
spatial and temporal correspondence not only with hard X-ray signatures, but also 
with enhancements in the white-light, leading researchers to also consider the role of 


pressure transients that are related to back 


chromospheric radiation (IDonea et al 


1999 


warming of the p 


lotosp here by the enhanced 


Donea fc Lindsev 


precipitation of electron (and p roton) beams fIZharkova fc Kosovichev 


more recently observed eve nts fj Alvarado-Gomez et al 


(jPedram fc Matthews 


2005). as wel 


as the direct 


20021) . However, 


20121 ) . and statistical studies 


20121 ) have indicated that the site of strongest white-light emission 


is not always coincident with the quake location. More recently, the role of the Lorentz 


force transients that occur whilst the coronal held is being reconhgured flHudson et al. 


20121 ). have been considered as another possible driver. Within this context one must 
also consider the fact that non -thermal particle beams carry a strong electric held 


fIZharkova &: Gordovskvv 


ambient plasma (IZharkova et al. 


20061 ). that in turn can induce an electromagnetic held in the 


20111 ) that modihes the magnetic held of the loop where 


particles precipitate. Magnetic held changes of this type are attractive as a potential driver 
since they do not require the same delay that the propagation of shocks requires, and there 
is increasing e v idence that many sunquakes are associated with erupting hux ropes (e.g. 


Zharkov et al. 


teOllal) : Green et al., (2015), in preparation). However, so far it has proved 


difficult to identify a dehnitive relationship between changes measured in the photospheric 
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magnetic field during the flare, and the site of the acoustic emission. In part, this may in 
fact be due to effects of the flare energy deposition, wh ich are known to produce artefacts in 


the magnetic held data 


here. 


Judge et ah 


flMartmez Oliveros et ah 


20141) . but in relation to the event studied 


fl2014j) estimate that the energy flux associated with the observed changes 


in the photospheric magnetic held are insufficient to account for the observed acoustic power 
of the sunquake. It is worthy of note that the sunquake in this event is somewhat unusual 
in that, while it is located in held that is of similar strength to the surrounding penumbra, 
there is no penumbra visible in this location in the continuum data. The majority of other 
quakes have been found in penumbral regions. 

Our investigations of the 29 March 2014 sunquake have focussed on spectroscopic 
signatures of plasma conditions in the overlying chromospheric, transition region and 
coronal layers, and their comparison with another location in the hare that di splays no 


acoust ic emission. These observations are complementary to those reported by 


Judge et al. 


(120141) which provide important insights into the photospheric environment around the 


sunquake, where we would expect to see signihcant energy deposition. Our holographic 
analysis indicates that the sunquake associated with this event is somewhat unusual 
in displaying a quite dispersed signal, but the spatial and temporal coincidence of the 
strongest egression kernels seen at 6 and 7 mHz and the time-distance source is compelling 
evidence for its reality. We hnd evidence for some variation in properties of the quake at 
diherent frequencies, with an earlier peak seen at 7 mHz than at 6 mHz. This hrst peak 
is coincident with the rise and peak of the flare impulsive phase, as observed in HXR 
between 17:45:31 - 17:47:46 UT (Fig. 1), while the second occurs during the decay of the 
HXR burst between 17:48:31 - 17:51:31 UT. As can be seen from Fig. 13, at 17:46 UT we 
observe a very close spatial coincidence between strong HXR emission at 25 - 50 and 50 
-100 keV, an intense, high velocity blue-shifted component in the hot Fe XXHI line (~280 
kms“^), strong red-shifted emission ( 70 kms in the Mg H 2791.6 A triplet line and 7 
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mHz acoustic emission. The He H 256.3 A line also shows down flows, but slightly offset 


from the other emission. The profile of the Mg H 2791.6 A trip 
core/wing intensity also appear consistent with simulations by 


et line, and t h e rati o of the 


Pereira et ah 


(1201511 that 


suggest a single strong temperature increase of 2000 K deep in the atmosphere. Taken 
together, these observations seem to strongly suggest the driver of the acoustic impact 
in the interval 17:47:46 - 17:50:01 UT is consistent with a thick-target scenario, where 
rapid heating by non-thermal electrons drives chromospheric evaporation and an opposite 
downward propagating shock that excites the sunquake. On the other hand, the later 6 
mHz peak in the interval ~17:51 UT seems to have a different driver, related to a second 
energy release at the SQ location at this time. In this case we observe no measurable HXR 
emission in the 25 - 50 and 50 - 100 keV range at the SQ location, but we observe very 
broad and red-shifted emission in all the Mg H lines (> 100 kms“^ in the k line). He H 
(~80 kms“Q, Si IV (~86 kms“^) and a small red-shift of ~26 kms“^ in hot (log T=7.1) Fe 
XXHI emission. The lack of HXR emission in the SQ location (this true at lower energies as 
well) is inconsistent in this case with the thick-target scenario, but the combination of very 
broad, red-shifted spectral profiles throughout the atmosphere may be consistent with the 
presence of a significant flux of Alfven waves. Using the measured line widths in the Mg H 
h and k lines we estimated the wave energy flux that would be associated with the presence 
of Alfven waves , and we compare this with the estimate of the peak acoustic power derived 


by 


Judge et ah 


(1201411 of 5 X 10® erg s“^ cm“^, both at 17:46 and 17:52 UT, finding values 
of 6.04x10^^ erg s“^ cm“^, and 3.87 xlO^^ erg s“^ cm“^ respectively. We stress that the 
impact of saturation at the earlier time, and our inability to actually measure the density, 
introduce significant uncertainties into these estimates, but we think it unlikely that these 
would be two orders of magnitude, suggesting that Alfven waves may indeed play a role in 
driving the sunquake in this event. 


We note that it is puzzling that we observe an apparently very strong chromospheric 
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response to the sunquake, while there appears to be a relatively weak response obser ved in 


the photosphere, particularly when the observations presented by 


Judge et ah 


fl2ni4h were 


so fortuitously located in space and time. We have no satisfactory answer to why this might 
be at this stage, but speculate that perhaps the wave energy flux that we observe (assuming 
for now that that is the source of the line broadening) dissipates between the regions of 
formation of the Mg II and Si I emission. What that might imply for the mechanism(s) 
of sunquake generation requires further investigation. Finally, we want to emphasise the 
complexity of the event studied here, which on the one hand seems to strongly support a role 
for chromospheric shocks in driving the acoustic disturbance, but also indicates the likely 
presence of a strong flux of Alfven waves. Coupled with the recent results that indicate an 
inconsistent correlation of sunquakes with the brightest white-light flare emission, and an 
association with the erupting flux rope, this highlights the difficulties in identifying a single 
driver, if indeed there is only one. 


The research leading to these results has received funding from the European 
Commission’s Seventh Framework Programme under the grant agreement No. 284461 
(eHEROES project). We would like to thank the referee, Phil Judge, for his thorough and 
insightful comments that have greatly improved the clarity of the paper. 
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Fig. 2.— Time-distance diagram for the acoustic source at (519 ,262 ) (left), with the 
theoretical travel time ridge over plotted (white dashed line) on the right. 
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Fig. 3.— IRIS Mg II slit-jaw image showing the location of the sunquake (SQ) and the 
non-quake (NQ) region and the location of the spectrograph slit (vertical line), with 6 (cyan) 
and 7 (magenta) mHz egression contours and RHESSI HXR contours (blue: 25 - 50 keV; 
black: 50 -100 keV). 
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Fig. 4.— EIS He II 256.3 A spectra at the location of the SQ (right column) and at the NQ 


location (left column) during the impulsive phase of the flare, as indicated by the vertical 
lines on the GOES and RHESSI HXR light-curves. 
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Fig. 5.— EIS He II 256.3 A spectra at the location of the SQ (right column) and at the NQ 
location (left column) during the gradual phase of the flare, as indicated by the vertical lines 
on the GOES and RHESSI HXR light-curves. 
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Fig. 6.— EIS Fe XXIII 263.7 A spectra at the location of the SQ (right column) and at 
the NQ location (left column) during the impulsive (black lines) and gradual phases (orange 
lines) of the flare, as indicated by the vertical ones on the GOES and RHESSI HXR light- 


curves. 
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Fig. 7.— IRIS Si IV A 1403 A spectra at the location of the SQ (right column) and at 
the NQ location (left column) during the impulsive phase of the flare, as indicated by the 
vertical lines on the GOES and RHESSI HXR light-curves. 
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Fig. 8.— IRIS Si IV A 1403 A spectra at the location of the SQ (right column) and at the 
NQ location (left column) during the gradual phase of the flare, as indicated by the vertical 


lines on the GOES and RHESSI HXR light-curves. 
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Fig. 9.— Top; Si IV intensities derived from spectral fitting for the SQ (left) and NQ (right) 
locations. Bottom: Doppler widths for the SQ (left) and NQ (right) locations 
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Fig. 10.— Mg II k and h spectra at the location of the SQ (right colnmn) and at the NQ 
location in the North flare ribbon (left colnmn) dnring the impnlsive phase of the flare, as 
indicated by the vertical lines on the GOES and RHESSI HXR light-cnrves. 
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Fig. 11.— Mg II k and h spectra at the location of the SQ (right colnmn) and at the NQ 
location in the North flare ribbon (left colnmn) dnring the gradnal phase of the flare, as 
indicated by the vertical lines on the GOES and RHESSI HXR light-cnrves. 
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Fig. 12.— Left: from top to bottom: Mg II k intensity; k/h ratio; Mg II k and h Doppler 
widths for the SQ location. Right: from top to bottom: Mg II k intensity; k/h ratio; Mg II 
k and h Doppler widths for the NQ location. 






-36- 


280 


270 - 


Hinode/EIS 29-Mar-2014 17:46:14.000 UT 


p 260 - 


250 - 


240 Lx 


7 mHz egression 
25-20 keV HXR 
50-100 keV HXR 
He II ri 
Mg II 2791.6 



500 510 

X (arcsecs) 


Hinode/EIS 29-Mar-2014 17:52:55.000 UT 



X (arcsecs) 


Fig. 13.— Left: Background image: Fe XXIII blue wing emission at 280 kms“^; 6 and 7 mHz 
egression contours in cyan/thick magenta; 25-50 and 50-100 keV HXR emission contours at 
80 and 90% of the maximum (purple and black, respectively); emission from He H at 60 and 
80 kms“^ (magenta), and emission from the red wing of the Mg H 2791.6 A triplet line at 
70 kms“^ (green). Right: as left, but at 17:52 UT. 
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